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Abstract 

In this talk I present the results from the measurement of the inclu- 
sive jet cross section and strong coupling constant based on the CDF 
Run IB data, and discuss prospects for Run 2. 

1 Introduction 

Jet production at hadron colliders provides an excellent opportunity for test- 
ing the theory of strong interactions, Quantum Chromodynamics (QCD). 
QCD has achieved remarkable success in describing hadron interactions 
at short distances (large momentum transfers), owing to the property of 
asymptotic freedom. Processes with large momentum transfer can then be 
described by an expansion in powers of a s (/•*)• Higher energy collisions 
produce jets of higher energy particles, while the value of a s decreases, 
improving the validity of perturbative expansion. Since the first measure- 
ments of the inclusive jet cross section by UA2, many experiments studying 
jet production published comparisons between theoretical predictions and 
experimental measurements. In particular, results from the CDF Run 1A 
data (1992-1993) inclusive differential jet cross section measurement, which 
showed significant excess of the data over the theoretical predictions at high 
Et, generated lot of excitement due to the possibilities of new physics, and 
motivated the re-evaluation of all components of theoretical calculations. 

The cornerstone of data to theory comparisons is the question of jet 
identification. In theoretical calculations, jets are manifestations of partons 
as relatively isolated sprays of energetic hadrons observed in the final state of 
high energy collisions. From the experimental point of view jets are defined 
as large energy deposits in a localized group of calorimeter cells, see Fig. [l]. 
To minimize the difference between parton level predictions and measured jet 

*Presented on behalf of CDF collaboration at " Corfu 2001: Summer School and 
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Figure 1: Lego display of jet event from CDF. 



properties, the clustering algorithms are used which could be implemented 
for both situations. Jet algorithms start from a list of "particles" that 
are taken to be calorimeter towers in experimental case and partons for 
theoretical calculations, and associate clusters of these particles into jets in 
such a way that the kinematic properties of the jets can be related to the 
properties of the energetic partons produced in the hard scattering process. 
This way jet algorithms allow us to see the partons in the final hadronic 
state. Differences in the properties of reconstructed jets when going from 
parton to hadronic or calorimeter level are of major concern for a good jet 
algorithm. 

Historically cone algorithms were used for jet measurements in hadron- 
hadron colliders. In Snowmass algorithms, two partons are combined into 
a jet if they fall within a cone of radius R in rj — <fi space, where R = 
\J {Arj 2 + A(p' 2 ) and Arj, A<j) are the separation of the partons in pseudo- 
rapidity and azimuthal angle. It is further required that the axis of the cone 
coincides with the jet direction as defined by Er weighted centroid of the 
particles within the cone, where Er = Esm.0 is the transverse energy. A 
similar algorithm with R = 0.7 is used in the experimental data analysis 
by using calorimeter towers instead of partons. The latest modification of 
this algorithm when applied to the theoretical calculations introduces an ad 
hoc parameter lZ sep || which is used to regulate the merging and separa- 
tion of parton clusters in manner similar to the experimental analysis. In 
addition to the difficulty with merging/splitting rules which were not ad- 
dressed in the original Snowmass algorithm, cone algorithms turned out to 
be sensitive to the soft radiation at the NNLO calculations. The second 
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Figure 2: (left) Variation in theory predictions for different renormalization 
scales; (right)The variation of the inclusive jet cross section for different 
TZ sep parameters. These calculations used the EKS program. 



type of algorithm, is so-called Kt algorithm which is infrared and collinear 
safe to all orders of calculations and gives close correspondence with jets 
reconstructed in calorimeter with jets reconstructed from partons, since it 
successively merges pairs of "particles" in order of increasing relative trans- 
verse momentum. The difficulties of applying this algorithm are due to the 
question of the subtraction of the underlying event energy and problem of 
minimizing computing time. 

There are additional issues which affect sensitivity of the measurements. 
First is the limited accuracy of the fixed order perturbative calculations due 
in part to the intrinsic uncertainties associated with the choice of renormal- 
ization and factorization scales, and fj,p, see Fig. §(left). For the usual 
range of fi = Et/2 to (i = 2Et the variation in prediction is around 20%. 

The variation in above mentioned lZ sep parameter attributes to some 
5-6% uncertainty, Fig. ^(right). But the limited knowledge or parton dis- 
tribution functions (PDFs), which are obtained from global fits to deep 
inelastic scattering (DIS), Drell-Yan production and other collider data, are 
the largest source of uncertainty for the calculations. The effect of this 
uncertainty will be discussed later. 
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2 Inclusive Jet Cross Section 



The inclusive jet cross section measurement is based on a data sample of 
integrated luminosity 87 pb~* collected by the Collider Detector at Fermi- 
lab (CDF) during the 1994-95 run (Run lb) of the Fermilab Tevatron pp 
collider operating at yfs = 1.8 TeV. The CDF detector is described else- 
where 0. Details of the measurement of the inclusive jet differential cross 
section can be found in Q. Briefly, jets are reconstructed using the itera- 
tive fixed cone algorithm mentioned above. The inclusive jet cross section 
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Figure 3: (left) Inclusive jet cross section from the Run IB data (1994 to 
1995) compared to a QCD prediction and to the published Run 1A data 
(1992 to 1993); (right) The ±lcr fractional change in cross section due to 
the dominate sources of systematic uncertainty. 

includes all jets in an event in the pseudorapidity range 0.1 < |??| < 0.7. 
The measured spectrum is corrected for the calorimeter response, resolution 
and the underlying event energy using an iterative unsmearing procedure 
which changes both the energy scale and the normalization simultaneously. 
Figure §(left) shows the corrected Run IB cross section compared to a QCD 
prediction and Run 1A results. As we see, the results of Run IB and Run 
1A are in good agreement. 

The majority of experimental systematic uncertainty arises from the un- 
certainty in the simulation of the response of the detector to jets. The 
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simulation is tuned to the data for charged hadron response, jet fragmen- 
tation , and 7T° response. Additional uncertainty is associated with the jet 
energy resolution, the definition of underlying event, the stability of the 
detector calibration over the long running periods and an overall normal- 
ization uncertainty from the luminosity determination. For each source of 
systematic uncertainty described above, except normalization, the inclusive 
jet cross section was re-evaluated by varying the corresponding parameter 
in the detector response by 1 a. For the normalization uncertainty it was 
changed by a scale factor, see Fig. |3] (right). 

Figure |I| shows the cross section compared to QCD predictions using 
three recent PDFs. Good agreement is observed using conventional (CTEQ4M 
and MRST) PDFs except at high values of transverse energy, starting at 200 
GeV, where an excess is observed. Similar deviation have been observed in 
other CDF jet cross section measurements such as the dijet mass ||, dijet 
differential [||, and T,Et || analyses. One possible explanation for the ex- 
cess of high Et jets is that the gluon distribution at high x is larger than 
conventional PDFs have indicated. We see that CTEQ4HJ, (special PDF 
set, which was created to explore the flexibility in the gluon distribution at 
high Et, and was generated by including CDF jet data in the global fit with 
higher statistical weight assigned to high Et points and a new parameter- 
ization of the gluon distribution ||), gives the best overall agreement with 
the data in overall shape and normalization. 



3 Strong Coupling Measurement 

The value of a s , a free parameter of QCD, is one of the fundamental con- 
stants of nature. Its determination is the essential measurement of QCD, 
and the observation of its evolution, or running, with momentum transfer is 
one of the key tests of the theory. 

The value of a s is determined by comparing the jet cross section with 
the next to leading order (NLO) perturbative QCD predictions [[?]]. In the 
Et region studied, the non-perturbative contributions to the inclusive jet 
cross section are estimated to be negligible ||. The procedure of extracting 
a s can be summarized by the equation 

^ = al(fi R )X^(fi F ,E T )[l + a s {fi R )k l (fx R ,^ F ,E T )] (1) 
where -42- is the transverse energy distribution of the inclusive jets, and 
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Figure 4: Run IB data compared to QCD predictions (EKS, /i = Er/2, 
■R sep = 1.3) using the CTEQ4M, CTEQ4HJ and MRST PDF's. Only sta- 
tistical uncertainties are shown on the data points. 



fip are the renormalization and factorization scales, a^(fiji)X^- \nF, Et) is 
the leading order (LO) prediction for the inclusive jet cross section, and 
a^ifi^X^fjip, E T )ki(n R , t*F, Et) the NLO contribution. Both X^ and 
h\ are calculated with the jetrad Monte Carlo program |J. All calcula- 
tions are performed in the modified minimal subtraction, MS, scheme. The 
jetrad Monte Carlo program has as an input parameters: jet cone radius 
R, TZ sep , PDF set and hr,/j,f- 

The inclusive jet data are divided into 33 Et bins, from which we obtain 
statistically independent measurements of a s for 33 different values of fiR. 
The a s values derived for (j,r = fip = Et using cteq4m || parton distri- 
bution functions (PDFs) are presented in Fig. || (left). For Et <250 GeV, 
there is good agreement with QCD predictions for the running of the cou- 
pling constant. The behavior of a s at higher Et values is a direct reflection 
of the excess observed in the inclusive jet cross section [^. 

The measured values of « s (//r) are evolved to the mass of the Z° boson, 
Mz, by using the solution to the 2-loop renormalization group equation. 
The values of a s (Mz) for all 33 measurements are shown in the inset of 
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Figure 5: (left) The strong coupling constant as a function of Et for hr = 
Et measured using cteq4m parton distributions. The shaded area shows 
the experimental systematic uncertainties. The curved line represents the 
NLO QCD prediction for the evolution of a s {Er) using a s (Mz) = 0.1178, 
the average value obtained in the region 40 < Et < 250 GeV. The a B {Mz) 
extracted from ol s {Et) is shown in the inset along with the weighted average 
as the horizontal line; (right) Uncertainties in a s (Mz) due to the renormal- 
ization scale [/,, (a), and parton distribution functions, (b). 



Fig. H(left). Averaging over the range 40-250 GeV, we obtain [ 10 1 

a s {M z ) = 0.1178 ± 0.0001(stat). 

Inclusion of the data with Et >250 GeV results in an increase of the average 
value by 0.0001. 

The experimental systematic uncertainties on the value of a s {Mz) are 
derived from those on the inclusive jet cross section and when summed in 
quadrature yield ^g^° uncertainty. 

The theoretical uncertainties are mainly due to the choice of renormal- 
ization and factorization scales and parton distribution functions. The scales 
Hp and hr are expected to be of the same order as the characteristic scale 
of the process, which in this case is the jet Et- We have evaluated the 
changes in a s (Mz) resulting from independently varying and fiR from 



Et/2 to 2Et and found that the largest changes occur for fiR = fj,p. For 
all results presented in this letter the two scales were set equal. The sen- 
sitivity of the measured a s (Mz) to changes in these scales is indicated by 
the shaded band in Fig. ||(right)(a). Over the Et range from 40 to 250 
GeV, the shift in a s (Mz) induced by changing the scale from Et/2 to 2Et 
is approximately ^4%, independently of Et- We use the cteq4a series to 
study the a s (Mz) dependence on the PDFs. In addition, we have studied 
a s (Mz) using PDF sets which do not include Tevatron jet results, such as 
the MRST(gt) set Q, the mrsa' series @, and two mrs-r sets JE3|. The 
X 2 , calculated by comparing the data with the theoretical prediction in the 
restricted range of 40-250 GeV, is used to quantify the agreement. The 
minimal \ 2 /d.o.f. = 1.38 is obtained for cteq4m (af DF =0.116), therefore 
we use this PDF in our final fit. Excluding the PDFs which have obvious 
disagreement {\ 2 /d.o.f . > 5), we estimate the uncertainty on the a s (Mz) 
due to the PDF choice to be ±5%, see Fig. fright) (b). 

The resulting value of a s (Mz) is in good agreement with the world 
average a s (Mz) = 0.1181 ± 0.0020. But what is more important running of 
a s is shown over a very large energy range in the single experiment. 



4 Dijet Triply Differential Cross Section 

The rapidity dependence of the dijet differential cross section provides infor- 
mation about new regions of x and Q 2 . Assuming a 2^2 hard scattering, the 
event kinematic variables (x,Q 2 ) are related to the jet's transverse energy 
Et and pseudorapidity r\ as 

X12 = Si^rexp^; Q 2 = 2E T 2 cosh 2 77* (1 - tanhr?*), 

where the sum is over all jets in the event. The x\,X2 are parton momentum 
fractions and Q 2 is the two body four momentum transfer of the interac- 
tion. CDF measured the dijet differential cross section [Q] where one jet 
is restricted to the central (0.1 <| r] |< 0.7) region and the other jet is 
constrained to different rapidity regions. Fig. ^ (left) shows the CDF data 
compared to the theoretical predictions. An excess similar to that of the 
inclusive jet cross section is observed. 
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Figure 6: (left)The dijet differential cross section for CDF from Run IB; 
(right) The inclusive jet cross section (number of events) in the central 
rapidity region for CDF. Predictions for Run lb, Run 2a, Run 2b using 
EKS NLO program and the CTEQ5M and CTEQ5HJ PDFs, as well as an 
extrapolation of the physics curves measured by CDF in Run lb. 



5 Run 2 Prospects 

One of the goals of the Run 2 is to obtain higher precision for QCD measure- 
ments. A data sample of 2/fo -1 will enable the jet cross section to be probed 
for higher Et and x values, than were possible at Run 1. In addition to the 
high statistics the increase in the center-of-mass energy from 1800 GeV to 
1960 GeV has a large effect on the high Et jet rate. For the full Run 2A 
sample the numbers of jets above 400 GeV are expected to increase from 11 
to ~ 500, with energy reach extending to 550 GeV, see Fig. ^ (right). For jet 
identification CDF will use both standardized jet cone and Kt algorithms. 

The Run 2 calorimeter upgrade provides calorimetry as precise in the 
forward region as in the central one. The inclusive jet cross section will be 
measured out to rapidity values of 3 in Run 2, and dijet differential mea- 
surement can be enhanced by allowing both jets to be non-central. In Run 2 
the measurement of the strong coupling constant will be possible from a Q 2 
of 100 GeV 2 to 300000 GeV 2 . With more luminosity delivered by Tevatron 
every day, we are looking forward to new QCD precision measurements at 
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CDF. 
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